Constraints on conductances for Y-junction of quantum wires 
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We consider the Y-junction, connecting three quantum wires, in the scattering states formalism. 
In the absence of fermionic interaction we analyze the restrictions on the values of conductance 
matrix, imposed by the unitarity of scattering S-matrix. Using the combination of numerical and 
analytical results, we describe the four-dimensional body of values of reduced conductance matrix. 
We show that this body touches the unit sphere at six points only, in accordance with Birkhoff-von 
Neumann theorem. It implies that the Abelian bosonization analysis for the vanishing interaction 
strength can be performed only in the vicinity of these six points. 



I. INTRODUCTION AND MODEL 

Investigation of the properties of electronic systems in 
reduced dimensionality is the subject of ongoing active 
studies, both theoretical and experimental. The effects 
of electronic interaction play a major role in one spa- 
tial dimension and lead to a formation of the Luttinger 
liquid stated The idealized Luttinger liquid is best de- 
scribed in terms of its excitations within the bosonization 
approach. The imperfections in quantum wires, such as 
impurities, are thoroughly investigated by different meth- 
ods, see^ and references therein. In view of potential 
technological applications, it is also desirable to theoret- 
ically analyze junctions of several quantum wires. Such 
analysis for the simplest variant of Y-junction connect- 
ing three wires was performed in a number of papers 
during the last decadeP^ These works divide into two 
groups. The first group considers the problem in the 
bosonization approach, which allows one t o fully take 
into account the bulk fermionic interaction! 3 ! 6 * 8 * ^ The 
shortcoming of this approach is the necessity to consider 
the Y-junction close to certain values of its transparency, 
as discussed below. Another group employs the con- 
ventional fermionic description for rather arbitrary Y- 
junctions! 4 * 5 * 12 * 13 ] and the drawback of these method is a 
finite order of fermionic interaction which can be taken 
into account. 

Recently, a method for the fermionic description of 
Luttinger liquid with a single impurity was developed, 2 
which performs a partial summation of fermionic dia- 
grams in the scattering states approach. This method 
allows one to obtain the exact scaling exponents known 
previously from the (Abelian) bosonization technique. It 
is possible to extend this method, formulated essentially 
in terms of non- Abelian bosonization, for the analysis of 
Y-junction} 1 ^ and to compare its predictions with those 
partly available from numerical and bosonization stud- 
ies. But before undertaking the study of the interact- 
ing case, it is necessary to clearly specify the restrictions 
on the transparency of Y-junction, which occur already 
at the level of free electrons and which are imposed by 



the basic quantum mechanical principles. The 3 x 3 S- 
matrix describing the Y-junction belongs to SU(3) group 
and consequences of it are sometimes not trivial. As we 
show below, the phase space for observable conductances 
(transparencies) is four dimensional (4D) and the val- 
ues fill the solid body of unusual shape inside the 4D 
unit sphere. Importantly, only six points of this body 
lie on the surface of such a sphere, which implies that 
the Abelian bosonization analysis can be done in close 
vicinities of these points only. 

This paper is organized as follows. We introduce our 
model and basic conventions in this section below. The 
relation between the S'-matrix and conductances is dis- 
cussed in Sec. [TT] The main results, numerical and ana- 
lytical, are presented in Sec. |III| We discuss our findings 



and make a comparison with previous studies in Sec. IV 



We consider a following setup: three semi-infinite non- 
interacting wires are connected by the Y-junction. The 
electron gas inside the wires is subject to transverse di- 
mensional quantization, so that only one channel per wire 
is available. We consider spinless fermions for simplicity. 
As usual, when interested in the response at smallest en- 
ergies, we can linearize the electronic dispersion around 
the Fermi level. In the present case, instead of defining 
right- and left-movers in each wire, it is more appropriate 
to speak of incoming and outgoing waves, with respect to 
Y-junction at the origin. The Hamiltonian of the system 
is given by 



H = vp 



^Lui* V ^W) (1) 



with the Fermi velocity vp set henceforth to unity. We 
unfold the setup in the usual wayJ^U^ by saying that 
the outgoing fermions on the negative £-axis are the 
fermions, which have passed through the junction, at pos- 
itive \x\, i.e. ip j>out (-\x\) ipj,out(\%\)- Thus, we have a 
three-fold multiplct of right-going (chiral) fermions. 

The boundary condition at the origin is described by 
the scattering S-matrix. For elastic scattering by the 
central dot, the outgoing fermions at the origin are con- 
nected to the incoming ones by the relation -0^ out — 
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In the scattered states representation, the right- and 
left- going fcrmionic densities acquire the form 

ipj,in(x) = tpj,out(x) = (S.ip)j(-x) 

^tin( x )^3,in(x) = pj(x), tp]^ (x)ljj^ out (x) = Pj{~x) 

(2) 

here and below we use the notation A = .A.S. 

For a multiplet of incoming fermions, VP = (ipi,ip2,'4 , 3), 
the incoming density pj = ^ -Pj-^S is given by diagonal 
matrix, 



(3) 

We also can write pi = \ (yf A o + ^|A 8 + A3) , p 2 = 

\ (\/i A ° + 73 A 8 - A s). h = TfAo - ^A 8 . Here the 
traceless Gell-Mann matrices, Xj, with j — 1,...8 are 
the generators of the SU(3) group, listed elsewhere (see, 
e.g.p^). In addition to these, we use also the matrix 

Ao = 1, proportional to the unit matrix, 1. We have 

the property Tr[\jXk] = 2Sjk, for j, k = 0, ... 8. 

The outgoing densities are given by pj = 
SjkSfaSkmiplipm = ^.S^.pj.S.^, so that in the 
matrix representation pj — S'.pj.S. We will mostly 
omit the hat sign over p below. 



II. S-MATRIX AND CONDUCTANCES 



A. Parametrization of the S-matrix 



that there are four independent components of the con- 
ductance matrix, and proceed with a reduced number of 
parameters in the description of S-matrix, Eq. ^ . 



B. Conductances 

The matrix of linear conductances Gjk connects the 
current Ij in wire j with the electric potentials \% in 
the leads through the relation Ij = J^kGjkVk- In lin- 
ear response theory^ the thermodynamically averaged 
currents are given by Ij(x) = j° oo dy((p j (x)p k {y)) - 
(pj(—x)pk(y)))Vk- In the static limit and in the absence 
of interaction the response functions may be evaluated 
to give 



Gjk = S jk - Tr(pjpk) = 5j k - \S jk f 



(6) 



One easily verifies that the charge is conserved, J^j Gjk — 
and that applying equal voltage to all wires produces 
no current, ^^Gjk = 0. 

In view of these conservation laws, it is more instruc- 
tive to discuss the current response to certain combina- 
tion of voltages. Let us define (I a , lb, Iq) = G.(V a , Vf,, V ), 
with 

V a = {V x - V 2 ), I a = (h - h)/2, 

V b = (V 1 + V 2 -2V 3 )/2 / b = (/ 1 +7 2 -2/ 3 )/3, (7) 

V = (Vi + V 2 + V 3 )/3, I = (h + I 2 + I 3 ), 

The connection between two matrices is given by G — 
A.Ftf.G.R.A, with A = diag{l/V2, ^/2~Ji, \/3) and 



The most general S-matrix is defined as follows 
S = I t 2 i, r 2 , t 2 3 

*31 ^32 r 3 



(4) 



where Tj is the reflection amplitude for wire j, and tjk is 
the transmission amplitude between wires j and k. The 
matrix S is unitary, S^S — 1, which allows its parame- 
terization via the exponent, S — exp ^ Ej=o ^3 A i) ■ ^ 
may be more convenient to use the Euler angles param- 
eterization of S. According td^, an arbitrary S'-matrix 
can be represented, up to an overall phase, as 



S = Ue l ^ e Ue tXaT , 

jj _ e iA 3 0/2 e iA 2 C/2 e iA 3 '0/2^ 



(5) 



with the range of most important parameters < £ < it, 
0<£<tt, 0<6»< tt/2, < ip < 4tt. 

Apparently, there is a redundancy in the description 
of S, as only the densities and not the fermion ampli- 
tudes enter the observable conductances. Below we show 



1/V2, 1/V6, 1/v^ 

u - I -1/V2, i/V5, i/V5 
0, -2/V6, i/Vsj 



(8) 



such that itfR = 1. 

It follows then that the third line and the third row, 
corresponding to Jo and Vq, are identically zero, Gjs = 
G^k = 0. We omit them for clarity below, and call the 
upper left 2x2 block of G by the same quantity G. These 
remaining four components are non-zero and we obtain 
the reduced conductance matrix in the general form 



G = 



G aa Gab 
Gba Gbb 



\ (l - iTr(A 3 A 3 ) 
-^Tr(A 8 A 3 ), 



-5j|Tr(AaAa) 
1 - iTr(A 8 A 8 ) 



(9) 



as is verified by direct calculation with Eqs. Q. 

For out purposes, it is more convenient to discuss not 
the reduced conductance matrix G, Eq.([9]), but the equiv- 
alent quantity M given by G = B.(l — M).B, with 
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B = diag(l/V2, 



Explicitly, we write 



M 



Tr(S^\ 3 SX 3 ) 
Tr(5tA 8 5A 3 ) 

z — x, y + T) 
y — T], z + x 



Tr(S*\ 3 S\ s ) 
Tr{Si\ a S\$) 



(10) 



We see that, due to the charge conservation, the ma- 
trix of conductances contains at most four independent 
components. Further reduction in the number of inde- 
pendent parameters occurs in case of additional sym- 
metries. For instance, the time reversal corresponds to 
G — > G\ so the unbroken time-reversal symmetry (T- 
symmetry) leads to G a b — Gba- The interchange of the 
wires 1 f> 2 leads to a change of sign for the off-diagonal 
components of G, so the symmetry between these wires 
leads to G a b — Gba — 0. The latter symmetry (1 O 2) 
in the presence of the magnetic field, however, leads to 
conclusion G a b = —Gba 7^ 0, since the parity change 
(123) — » (213) should be accompanied by the change of 
sign of the magnetic flux piercing the Y-j unction, and 
thus to transposition of G. 

These arguments are obviously applicable to a junction 
of n wires, i.e. the number of relevant physical parame- 
ters is (n — l) 2 , which is the dimension of the reduced 
conductance matrix. T-symmetry leads to symmetrical 
form & — G, thus reducing the number of parameters 
to n(n- l)/2. 



C. Reduced parametrization 

The general expression ^ contains eight Euler an- 
gles; however, the observable conductances ^ include 
only traces of products A 3 ( 8 ) an d S'X 3 ^S. Similarly, 
in the presence of fermionic interaction, each correc- 
tion to conductances contains only products of A 3 ( 8 ) and 
SU 3(8) SP Since A 3 and A 8 commute with each other, 
the angles <fi, -ip,T in Eq. ^ drop out of any observable 
quantity. Further, the explicit calculation shows, that 
only the combination if) + cf> enters our formulas and we 
can redefine ip — > t/j — cf> in ^ (or, equivalently, setting 
<t> = above), to eliminate (f> from the subsequent analy- 
sis. (More precisely, cf> may be generated by the RG flow, 
but it does not enter the right-hand side of RG equations 
and is absent in the conductances.) Hence, without a loss 
in generality, we may parametrize 

S _ e iA 2 f/2 e iA 3 (7r-^)/2 e iA5<yA 2 e72^ Q-n 

For the above matrix M we have 

M11 = \ cos £ cos £(1 + cos 2 9) 4- sin £ sin £ cos 9 cos ip, 



M12 = cos £ sin 2 6>, M 2 i 
M 22 = |(l + 3cos20). 



= — ^ cos £ sin 2 ( 



(12) 

From Eq. ([12) it is clear that -1/2 < M 22 < 1 for arbi- 
trary S'-matrix. It follows from Eq. ([9| that the tunneling 
conductance < Gbb < 1, as it should be. 



In the next section we consider two important partic- 
ular cases. First is the time-reversal symmetrical case, 
when we should have rj = in (10). Another case, partly 



considered irPHSIISI corresponds to symmetry between 
the wires 1 and 2 and the magnetic flux piercing the Y- 
junction, which reads as y = in Eq. (10). 



It should be noticed, that the correspondence between 
the arbitrary S-matrix of the form ( |11[ ) and the under- 
lying microscopic Hamiltonian is not simple beyond the 
lowest-order Born approximatioiM It means that the 
effective low-energy description for a particular micro- 
scopic setup can produce further restrictions on S, which 
add to those stemming from the above general symme- 
tries. 



III. ALLOWED VALUES OF CONDUCTANCE 

A. T-symmetric arbitrary Y-junction 

We showed that the values of conductance are uniquely 
defined by the matrix M. Let us discuss the case rj = in 
(10), which co rresponds to time-reversal symmetry and 
to £ = £ in ( 12 1. Let us take an arbitrary S'-matrix, given 
by a random realization of ip,9,^ in (12), and for each 



matrix numerically calculate three numbers (x,y,z) as 
defined by ( 10 ). The possible values (x, y, z) do not fill a 



cube or a sphere, but rather form a peculiar tetrahedron, 
which is shown in Fig. [T] 

The projection of this body onto the (x, y) plane is 
shown in Fig. [2^,. The limiting curve, bounding the 
possible values (2, y), is called deltoid curve and is given 
parametrically by 

x = (2 cos t + cos 2t) /3, y = (2 sin t — sin 2i) /3, 

(13) 

with t e (0,2tt). 

Let us prove, that the shape in Fig. [2^, is a deltoid. 
Introduce the notation r — cos 9, u — cos£. Clearly, for 
fixed r, u the boundary is reached at cos ip = ±1. In view 
of the symmetry (tJj — > n — ip, 9 —¥ 7r — 9), we may choose 
cos ip = — 1 and seek the extrema of M' for both signs of 
r. We have 



x = -|(1 + r)(l + u 2 + t(-3 + u 2 )), 



y = ^(-l+r>. 



(14) 



We find the extremum curve by i) fixing the direction 
d{x jy) — which gives the condition between du and dr 
and ii) expressing du via dr from this condition, we find 
the maximum of x, (simultaneously of y). It produces 
the parametric description of the limiting curve 



x=i(3r 2 -l), 2 / = ±i(-l + r) v / 3 + 6r-9r 2 , 

with the range of r g (— 1/3, 1), defined by positiveness 
of the square root in ( |15[ ). This curve coincides with ( 13 ) 
after substitution t = (1 + 2cost)/3. Importantly, the 



FIG. 1: (Color online) Three views of the three-dimensional 
body, representing all possible values of the "conductance" 
matrix M in the situation with the unbroken time-reversal 
symmetry. 



above relation x — |(3t 2 — 1) means that z = (i.e. 
Mn = M22) for the limiting curve. 

The cross section y — is shown in Fig[2]3. The limit- 
ing curves are straight line x + z = 1 (corresponding to 
zero conductance Gbt), and a parabola, given by 

x = (3i-l)(* + l)/4, z = (3i + l)(t- l)/4. (16) 

The lowest point x = y = 0, z = — 1/3 corresponds to 
maximum transparency of K-junction, discussed below. 

The projection onto the (x, z) plane is shown in FigjlJ;. 
This projection in fact corresponds to two cross-sections, 
delimited by parabola and straight line, and rotated by 
an angle 27r/3 in three-dimensional space. The resulting 



FIG. 2: (Color online) (a) A top projection of the "tetrahe- 
dron", showing the deltoid curve, denoted by a red line. The 
circumscribed circle is also shown, (b) The section y — 0, cor- 
responding to the symmetry between first and second wire; 
the bounding curve consist of a parabola and a straight line, 
(c) The projection onto the (a;, z) plane; one can see parabolic 
and straight lines, bounding the body, similarly to the cross 
section y = 0, and shown by red lines. 



contraction by a factor of sin(27r/3) = \/2>/2 along the x- 
axis is taken into account, when drawing limiting curves 
here. 

Notice that there are only four points of "tetrahedron" 
which lie on the surface of the unit sphere in Fig. [Tl They 
are N = (0, 0, 1), A 3 = (1, 0, 0), A 1>2 = (-1/2, ±73/2, 0) 
in coordinates (x, y, z) and correspond to three detached 
wires, N, and (three times by) one detached wire, Aj, 
respectively. The related matrix IS^fcl 2 in ^ in cases 
N, Ai, A 2 , A3 has the form of permutation matrix, Eq. 



( 20 1 below, with unity value for po , Pi , P2 , P3 , respectively. 
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B. Magnetic flux and symmetry between two wires 

In case of magnetic flux piercing the Y-junction and 
the additional symmetry between the wires 1 and 2, we 
have the relation y — in in Eq. ( 10 ). At the formal level, 



however, this case is fully analogous to the one considered 
above. Indeed, the previous case corresponded to £ = £ 
and the present case reads £ — ir — £. This difference 
amounts to a change (x,y, z) — > (z,rj,x), which means 
that the Fig. [I] and Fig. [2] having the same form in new 
axes. 

The equations (14 1, (|15[) are now valid for z, rj, and 



the limiting curve, bounding the values of conductance 
matrix, is again deltoid. Further, it happens at x = 
0, i.e. for the totally symmetric situation with respect 
to permutation of all three wires. Such particular case 
of total symmetry was considered earlier for interacting 



fermions irPE 

Similarly, there are four points which lie on the sur- 
face of the unit sphere in coordinates (z, r\, x). They 
are the previous points N, A3 and new ones x± = 
(-1/2, ±^3/2,0). The latter points are the chiral fixed 
points corresponding to currents flowing from the wire j 
to j =F 1 and discussed irP. 



C. Geometrical interpretation 

In this subsection we discuss the restrictions for the 
elements Mjk, which stem solely from the unitarity of S- 
matrix and do not involve particular properties of SU(3) 
group. We represent Eq. (|17[) in the form 



M 



(A 3 |A 3 ), (A 3 |A 8 ) 
(A 8 |A 3 ), (A 8 |A 8 ) 



a cos a, 6cos/3\ 
a sin a, 6sin/3 J ' 



. (17) 

here the shorthand notation (x\y) — ^Tr(x.y) was intro- 
duced. 

It is seen, that the columns of M are projections of A3 
and A 8 onto the hyperplane spanned by A3 and A 8 . We 
can think of the vectors a — A 3 (A 3 |A 3 ) + A 8 (A 8 |A 3 ) and 
b = A 3 (A 3 |A 8 )+A 8 (A 8 |A 8 ). Each of the sets {A,} and {A^} 
with j = 0, . . . 8 forms the orthonormal basis for Ilerrm- 
tian matrices with the scalar product (-|-). Therefore we 
have a 2 = |(A 3 |A,)| 2 < £-=i |(A 3 |A,)| 2 = 1 and 

similarly b 2 < 1. Notice that these inequalities follow 
from the unitarity of S, i.e. charge conservation. 

Further, let us consider a unit vector v = cos7A 3 + 
sin7A 8 in the hyperplane spanned by A 3 and A 8 . The 
projection of the circle 7 G (0, 2tt) onto the plane {A 3 , A 8 } 
is given by «p = cos 7a + sin 7 b and should have the 
form of the ellipse. A simple calculation shows that the 
principal axes of this ellipse are given by 



cos 2 C± = \{a 2 + b 2 )±\^{a 2 - 6 2 ) 2 + 4(ai) 2 , (18) 

here a.b = abcos(a — /?) Evidently, £_ is the dihedral 
angle between the hyperplanes {A 3 , A 8 } and {A 3 , A 8 }. We 



also have an obvious geometrical inequality, cos 2 £+ < 1 , 
which can be written in the form 



a 2 + b 2 < I + a 2 b 2 sin 2 {a- (3) 

s to the statement thi 
the non-negative matrix M^.M = 



(19) 



and amounts to the statement that the eigenvalues of 

a 2 , a.b^ 



a.b, b 2 



, are less 



or equal to unity, ||MtM|| < 1. 



D. Birkhoff-von Neumann theorem 

Introducing Mjj = |SV/| 2 in Eq. (jfjj), we write Gij — 
Sij - Mij. We have the property £ l M ij = 12 j Mij = 1, 
the matrices of this form are called doubly stochastic. 
For a vector space with elements X t we may form a lin- 
ear combination of the form J2i c i-^i- If each Ci > 
and Ci — 1, then such combination is called convex 
combination. 

The Birkhoff-von Neumann theorem states that an n x 
n matrix over 3? is doubly stochastic if and only if it is a 
convex combination of permutation matrices. 

A doubly stochastic matrix is called unistochastic if 
its elements can be represented as squared moduli of el- 
ements of unitary matrix. Any unistochastic matrix is 
contained in the set of doubly stochastic matrices, but 
not vice versa. In our case a general doubly stochastic 
matrix is a candidate for conductance matrix, but not 
all such matrices have a quantum-mechanical counter- 
part, S— matrix, i.e. they do not belong to unistochas- 
tic set. The analysis of 3 x 3 unistochastic matrices, 
done in [20 21 , is relevant to our study. Particularly, 
it was shown, that the extreme points of unistochas- 
tic set are permutation matrices. The doubly stochastic 
property represents the classical Kirchhoff's rule for the 
charge conservation, whose quantum-mechanical coun- 
terpart, the unitarity of S— matrix, necessarily involves 
the complex-valued quantities, except for extreme points 
and surfaces. A special role is played by two subsets of 
matrices of the form (cf. Eq.(16) hP^}) 



'PO, P4, P5,\ I Pi, P3, P2,' 

Ml = I Pb, PO, Pi, I , M 2 = p 3 , P2, Pi, 
\Pi, Pb, PO I \P2, Pi, P3 

(20) 

The above mentioned extreme cases correspond to pi = 1 
for a given i, with all other pi = 0. These cases corre- 
spond to the above points N, A\, A 2 , A 3 , \+, X-, f° r 
i = 0, . . . 5, respectively. 



Explicitly, we have in general case of M, Eq. ( 10 ) 

x = P3 - +P2), y = ^-(pi-ps), 



z=PQ-\{pi+Pb), r)= 2~(P±-P5) 



(21) 



with Y,i=oP* = L 



6 



The Birkhoff-von Neumann theorem does not imply a 
uniqueness of the convex combination. This becomes ev- 
ident when we notice that Mi + M 2 = for arbitrary p , 
if Po = P4 = P5 = —pi = —P2 = —P3 ■ For further illus- 
tration of this issue, consider the above additional con- 
dition of unbroken time-reversal symmetry 77 = which 
reduces the freedom in choice of parameters Pj , but not 
completely. The allowed values of x,y,z for the doubly 
stochastic matrix Mi + M2 of the above form fill the 
triangular bipyramid, which is the "classical" analog of 
the body shown in Fig. [T] This bipyramid is defined by 
its vertices (x,y,z) at (0,0,1), (1,0,0), (-l/2,±V3/2,0) 
and (0, 0,-1/2). The consideration of the case with mag- 
netic flux, y = 0, 77 7^ is done similarly to above, due to 
obvious symmetry of Eq. (21 1. We have the same bipyra- 
mid in coordinates (z, 77, x), as depicted in Fig. [3] 

Concluding this section, we make two observations. 
First, there are only six points, which lie on the surface of 
the unit four-dimensional sphere (x, y, z, rj). They are N, 
^1,2,3) X± an d this set is the same for doubly stochastic 
and unistochastic matrices. Second, the above inequality 
( 19 1 can be cast into the form 



yV + y 2 + \/z 2 + if < 1. 



(22) 



This condition at r\ = represents two cones. Three 
straight lines, connecting N and Aj lie on the upper cone 
z = 1 — \J x 1 + y 2 and belong to the boundary of allowed 
values of conductance both for doubly stochastic and 
unistochastic sets. The lower cone, z = — 1 + y 1 x 2 + y 2 7 
does not touch the above "tetrahedron" in both variants. 



IV. DISCUSSION 

A. "Stability of the action" 

Consider the energy, associated with the chiral 
fermionic densities in wires. For the ring of length I the 
energy of the chiral current pj around the energy min- 



imum is given by irvpl 1 (pj 



off, see e.gPUl and 

references therein. Here p^ is discrete variable (zero 
mode), measuring the number of particles moving at a 
given direction. Extending this formula to our situation, 
we have the energy of the incoming and outgoing d.c. 
currents per unit length 



E in = ttvf 



3 ,0 

2J (r 1 / dx(p jM (x)) 

A 1 J —I 



3 



E ut = nv F y2 (l 1 dx {pj,out(x)) 



(23) 



with I — > 00. This formula is written for the non- 
interacting fermions, and with some restrictions it can 
be used in the interacting CclSG, cLS well. Such generaliza- 
tion is possible, e.g., when one considers the finite region 




- 1 




-05 oa 




FIG. 3: Three views of bipyramid, representing all possible 
values of classical "conductance" matrix M in the symmet- 
ric situation (y — in Eq. (21 1) with broken time-reversal 
symmetry. 



of the interaction close to the origin, |x| < L, in Eq. 
( 23 1 . In this case the contribution of the interacting re- 
gion x G (-L, L) in the semi-infinite leads to Ei ntOUt is 



negligible, and we can adopt the single-particle descrip- 
tion of the outgoing states and the free-fermion relation 
between the applied voltage and the incoming current, 



(P:j 



-Vj. The existence of the non-interacting 



leads beyond the finite region of interaction is impor- 
tant for the proper definition of S'-matrix; it also leads to 
the absence of vertex corrections to conductance and the 
(correct) unity value of conductance of clean Luttinger 
liquid, as discussed at length iiP. Here (pj.in) stands 
for zero mode of the current, i.e. (pj,in(x)) averaged over 
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x in (23). The linear response theory gives for the net 
current Ij = (pj,i„) - (pj,out) = \G jk V k = G jk (p kiin ). 



Therefore (pj jOUt ) = M jk (p kiin ), 

Mj k = Sjk — Gjk = |5*jfc| 2 - Notice that, by construction, 
Tr(M) > for any S. 

The energies ( 23 ) represent only a part of the energy 



fermion, ipj( x ) ~ tpj(x) expikpjX, with Fermi momenta 
kpj = kp + Vj/vp; here we extracted the "smooth" part 
with the above matrix °f ^ ne f errmon operator 'tpj(x), see e.gPS. We see now 



of chiral fermions associated with the zero modes of the 
system. Indeed, from the general inequality (J fg) 2 < 

(I P)(1 9 2 ) we obtain irv F fl' 1 f° ; dx (pj, ln { x ))) ^ 

Trvpl^ 1 J_ ; dx (p 2 in (x)} = E to t, where the last quantity 
is the total energy density in bosonization description. It 
is assumed, that the total energy of incoming fermions is 
given by the zero modes, EV„ = E to t- 

Since the scattering is elastic, and a part of the energy 
Ei n after the scattering may become distributed among 
the excited particle-hole pairs (with non-zero momentum 
in appropriate description), one obtains the inequality 
Ei n _ E out which was called a stability condition irM 
This inequality means, that the eigenvalues of the non- 
negative matrix M*M, or equivalently M^M, are less or 
equal to unity. 

We saw above, however, that the stability condition 
||MtM|| < 1 follows automatical ly from the unitarity 
of S. The bosonization approac h 1 6 ! 7 ! 18 * 19 ! does not as- 
sume the unitarity of S'-matrix (see below), but rather 
starts from the consideration of conductance matrix, G, 
or "current splitting" matrix, M . As shown above, Eqs. 
(17), (19), the restrictions for the parameters are 



a 2 < 1, 



b 2 < 1, 



cos 2 (a -0)< (a- 2 - l)(b- 2 - 1) 



(24) 



The case considered in Ref.! 6 -! reads in our notation as 
1 - |G S , asina = |A, b = a, (3 = a + n/2. We 



a cos a 



see that the unitarity-related inequalities (24) reduce in 



this case to (1 — §G S ) 2 + (§ A) 2 < 1, which correspond to 
Eq. (11.13) irP. At the same time, Figs. [2] [3] show that 
the bounding estimates p4|), (22) are too weak. 



B. Abelian bosonization 

Let us further comment on the nature of the inequal- 
ity E in > E out linking it to the bosonization approach. 
According to scattering states formalism described above 
and in Ref. [2 the components of outgoing density may 
include non-diagonal interference terms. For example, 
the appearance of A4 component in the outgoing cur- 
rent pi indicates the presence of combination ifypx+h.c. 
Contrary to the interpretation irfSl, the existence of such 
components of the density does not mean the power dis- 
sipation, as the scattering is elastic. The difference be- 
tween the diagonal terms ipi^j and off-diagonal terms 

ipjipi (i 7^ j) becomes more transparent, if we consider a 
non-equilibrium situation with voltage biases Vj 7^ 0. In 
this case we should recall the existence of the rapidly os- 
cillating exponents in the wave-function of the incoming 



that the diagonal terms ipj^j m t ne outgoing current 
contain only smooth part of the density, whereas the off- 
diagonal part of the density contains "beats" in its pro- 
file, •tjy^i expi(Vj - Vi)x/v F , cf. RefPSl. Evidently, the 



off-diagonal terms vanish upon the averaging in ( 23 ) , and 
this is exactly when the energy difference Ei n — E out oc- 
curs. This difference is thus associated with the weight of 
the off-diagonal components of fermionic density, which 
may also include an analog of 2k F Friedel oscillations. 

The condition of exact equality Ei n = E out indicates 
a possibility to fully express the outgoing currents in 
"smooth" components of incoming chiral densities. It 
means that, in the appropriate basis, the whole Hamilto- 
nian is expressed entirely in smooth densities which is the 
essence of Abelian bosonization approach. Notice, that 
it is always possible to express the kinetic part of the 
action in terms of smooth densities in each of the semi- 
infinite wires. The difficulty in such a description arises 
at the moment of imposing a boundary conditions, which 
connect the smooth densities in different semi- wires ]SH2i 
Thus, for instance, a simpler case of the wire with a sin- 
gle impurity allowed a simple perturbative analysis only 
in two limiting cases of a clean wire and two detached 
semi- wires Psl 

Evidently, the condition Ei n = E out happens at a = 
6=1 and a = j3 ± it/2 in (17). We have two classes 
discussed ir i 8 l n l 27 l namely 



Mi = 

M 2 = 



cos a, sin a 
-sin a, cos a 



- cos a, 
sin a, 



sma 
cos a 



(25) 



(the form of M\ in our notation has a direct corre- 
spondence with Eq. (10.17) irP ). In these cases M 
is orthogonal, M'M ~ 1, and one may promote the 
above equality for the averages to the operator form, 
Pj,out = MjkPk,in- Then the commutation relations for 
the chiral fields, \fij,i n (x), p k ,in{v)] = ^S jk S'(x-y), pre- 
scribed in bosonization, are fulfilled also for Pj y0 ut- 

In (Abelian) bosonization approach, the consistent 
perturbative analysis can be done only in a close vicin- 
ity of the RG fixed point, which is scale invariant or 
conformally invariant point in iD.EEffl The require- 
ment for the above matrices ( [25] ) to represent the RG 
fixed point in bosonization leads to additional boundary 
conditions for three chiral densities at the Y-junction. 
These boundary conditions are expressed via a set of 
Neumann and Dirichlet boundary conditions for appro- 
priate bosonic fields. An important observation, made 
in Ref. [T5], is that the trace of the matrix M at the 
fixed point is an integer number and equal to the num- 
ber of Neumann conditions imposed. In simple words, 
Tr(M) = 1 + Tr{M) = 1 + 2z, at the fixed point should 
be equal to the number of wires, detached from the Y 
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junction; it becomes rather obvious when we recall that 



Tr(M) = M 2 + H 



r 3 \ 



Evidently, the number of 



Based on the above observations I make a following 
stronger conjecture. The two classes (25) do not rep- 



detached wires can be 0, 1,3, but not 2. 

The fixed points discussed in bosonization approach in- 
clude the above N point, Mi (a = 0) = ( J' x ) , for three 

detached wires ; the A 3 point, M 2 (a = 0) = (~ ' i), 
for the third wire detached, and also Mi(a — ±27r/3) = 
( , ±v ^/ 2 ) for situations with the first and second 

\±V3/2, —1/2 J 

wires detached, points A\ and A2. Two chiral fixed 
points \±, present in the broken time- reversal case, are 

M 1 (a = ±2tt/3) = (~$ 3> ^f/f); we have zero de " 
tached wires in this case. 

However, formally there is seventh fixed point, first 
discussed irP and denoted as Dp irPtSLLU It is given 
by Mi (a = 71") = ( V' From our Eq. we have 

G aa — 1, Gt,b — 4/3 > 1, i.e. the value of one-terminal 
conductance Gbb at the Dp point is larger than the 
maximum value of unity for the non-interacting situa- 
tion. This enhanced value of conductance was interpreted 
as the signature of Andreev reflection and multiparticle 
scattering at the Y-junctionPSl We observe that for Dp 
point one has Tr(M) = —1, and it contradicts the charge 
conservation law both in quantum mechanical and clas- 
sical formulation, at least for free fermions. Indeed, we 
should have Tr(M) > for any unitary S-matrix, and for 
any doubly stochastic matrix, representing Kirchhoff's 
rule. 

One might argue, that our analysis, performed for free 
fermions, is inapplicable to the case of strongly interact- 
ing regime, where the impor tan ce of Dp point was em- 
phasized. It was predicted irP^that the fixed point Dp 
becomes stable at strong attraction between fermions, 
g > 3, whereas it is unstable at weaker interaction, g < 3 
(cf. also^). The existence of Dp point, however, was 
not ruled out for any values of interaction and, particu- 
larly, in the free fermion limit, g 1. The inspection 
of formulas (10.29)-(10.31) irP and Eq. (62) irP for the 
scaling dimension of leading perturbations confirms it. 
Therefore we have to conclude, that the discussed bound- 
ary for physical conductances was not so far detected 
by means of conventional bosonization. One possibility 
here, potentially favorable for the existence of Dp point, 
is that this boundary changes with the Luttinger param- 
eter g and thus should be separately discussed for each 
v alue o f g. Another possibility, supported by our analysis 
in™, is that the boundary determined for free fermions 
remains the same for any interaction strength. Particu- 
larly, any physical RG fixed point belongs to this bound- 
ary, while the above Fig.^ and the Fig. 5 irP^l depict the 
same parabola. From this latter viewpoint, the discussed 
Dp point lies outside the physical domain and is not con- 
nected with the inside of the body in Fig. [T] by any RG 
flow. As was already mentioned in Sec. |IVX| we assume 
the infinite non-interacting leads beyond the finite region 
of interaction, which setup provides the proper definition 
of S'-matrix and vanishing vertex corrections to conduc- 
tance. 



resent any physical realization of Y-j unction, except for 
six points listed above and corresponding to permuta- 
tion matrices in Eq. ( 20 ) . These six points are the ex- 



treme points of the four-dimensional body, whose three- 
dimensional projections are shown in Fig. [T] All other 
physical realizations of the Y-junction lead to values 
{x,y,z,rj) lying strictly inside the four-dimensional unit 
sphere, Eqs. (22 1, ( 24 ) , whose surface is needed as a start- 



ing point in bosonization. Thus the Abelian bosonization 
study of the Y-junction can be performed at six RG fixed 
points only, at the vertices of the body in Fig. [I] If the 
RG fixed point occurs at the boundary of this body in- 
stead, then it is unsuitable for conventional bosonization 
analysis. It should be stressed again, that this conjecture 
is based on the present free-fermion analysis and on the 
results of perturbative two-loop RG treatment valid for 
relatively weak interaction.^ The case of strong interac- 
tion is to be discussed elsewhere.^ 



C. Larger-TV junctions 

The observations reported in this paper for the Y- 
j unction, can be generalized for a larger number of semi- 
wires, N > 3, attached to one spot. The S'-matrix in this 
case will be given by a finite rotation matrix in SU(N) 
group, and the squares of its entries form the unistochas- 
tic matrix M. The set of all such M lies within a set of 
bistochastic matrices, whose extreme points are N\ per- 
mutation matrices, according to the Birkhoff-von Neu- 
mann theorem. These permutation matrices describe the 
physical situations when some of the wires are detached 
from the central spot, whereas the other semi-wires form 
either ideal wires (in case of pairwise partial permuta- 
tions) or chiral groups, involving three or larger n umber 
of semi-wires. For example, for N = 4 (see also 1 20 * 21 *) 
we obtain 24 matrices, which correspond to 4 detached 
wires (one point in phase space), 2 detached wires and 
one full wire (six points), one detached wire and a chi- 
ral group of three semi- wires (eight points), two full wires 
(three points), and a chiral group of four semi- wires (3!=6 
points). These points lie on an appropriate unit hyper- 
sphere and thus are suitable as fixed points in Abelian 
bosonization. We note also that the particular case of 
fully symmetric S-matrix and fully symmetric interact- 
ing wires can be solved within TBA approach^! for any 
N, because the high symmetry effectively reduces this 
case to N = 2. The latter approach allows to consider 
the RG fixed point, corresponding to the perfectly trans- 
mitting junction [z = — 1/3 in our above case N = 3). 
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